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1 Introduction 

The decomposition of the nucleon spin projection of 1/2 (in units of %) into contributions from helicities 
and orbital angular momenta of partons became a topic of major interest in experimental and theoretical 
hadron physics after the European Muon Collaboration at CERN had published the surprising result 
that quark helicities contribute only an unexpectedly small fraction [T]. Extensive nucleon spin studies 
were carried out at CERN EH, SLAC H, DESY Q, JLAB © and at BNL 00. From the data, 
parton helicity distributions of the nucleon were extracted using the framework of perturbative QCD. By 
now, the contribution of the quark helicities to the nucleon spin is known to be about 30%, significantly 
smaller than the value of 60% expected from the Ellis-Jaffe sum rule [9]. Relativistic quark motion is 
responsible for the reduction from the value of 100%, expected in the naive quark-parton model iPTOl . In 
spite of the ongoing theoretical debate on how to correctly perform a gauge-invariant decomposition of 
the nucleon spin, agreement exists that besides the contributions of the quark helicities also the gluon 
helicity contribution AG is a measurable, gauge-invariant observable (see, e.g. Ref. [1 1J). 

The data of present inclusive polarised deep-inelastic scattering experiments cover a narrower range in 
the photon virtuality, Q 2 , as compared to unpolarised ones and hence their QCD analyses (e.g. Ref. [12]) 
show only limited sensitivity to the gluon helicity distribution as a function of the gluon momentum 
fraction Ag(x), and to its first moment, AG. Such a determination of Ag(x) from QCD evolution has 
therefore to be complemented by direct, dedicated measurements. 

Direct determinations of the average gluon polarisation in a limited interval of x, {Ag / g) , were performed 
in a model-dependent way using the Photon-Gluon Fusion (PGF) process by SMC Ifl3l , HERMES Ifl4l 
and COMPASS 03]]. These analyses used events containing hadrons or hadron pairs with high transverse 
momenta, typically 1 to 2 GeV/c. This method provides good statistical precision but relies on Monte 
Carlo generators simulating QCD processes. PYTHIA 021 was used by HERMES and by COMPASS 
for the analysis of small Q 2 events, and LEPTO H3 by SMC and COMPASS for Q 2 > 1 (GeV/c) 2 
events IPT81 . All measurements yield a small value of the gluon polarisation at x rj 0.1. This is consistent 
with recent results from PHENIX [7] and STAR [8] at RHIC, where the production of inclusive 7T° or 
high transverse momentum jets led to constrain the magnitude of (Ag/g). 

In this paper, we present new results on {Ag/g} and the virtual photon-nucleon asymmetries obtained 
from charm production tagged by D meson decays in 160 GeV/c polarised muon-nucleon scattering. 
The data were collected by the COMPASS Collaboration at CERN in the 2002-2004 and 2006-2007 
running periods. The results supersede the ones given in Ref. [ 19'] since they are based on the full data 
sample and an improved analysis method; additional final state channels are added as well. The gluon 
polarisation is determined assuming that open charm production is dominated by the PGF mechanism, 
y*g — > cc, as depicted in Fig. [T] The subsequent fragmentation of the cc pair, mainly into D mesons, is 
assumed to be spin-independent. The dominance of the PGF mechanism in the COMPASS kinematic 
region is supported by the EMC results on F 2 CC (Ref. |[20l . further discussed in Ref. EH), and by a 
COMPASS study of charm meson production 1221 . The determination of the gluon polarisation based 
on this assumption, although limited statistically, has the advantage that in lowest order of the strong 
coupling constant there are no other contributions to the cross section. 

In the present analysis, only one charmed meson is required in every event. This meson is selected 
through its decay in one of the following channels: D*(2010) + — > D ?!^ — > (K _ 7r + /K _ 7r + 7r°/K _ 7r + 
n + n~)n^ ov/ or D — > K~n + , as well as their charge conjugates. The former samples are called 'tagged' 
ones while the latter is denoted 'untagged'. Virtual photon cross section asymmetries, ^r N ^ D ° x ; an d the 
average gluon polarisation {Ag/g) are extracted from these open charm events. In Table[Tj the kinematic 
variables describing the /J.N scattering process are listed. In this analysis we have also employed, for 
the first time, next-to-leading order QCD calculations for the determination of the gluon polarisation. 



Throughout this paper, x denotes the gluon momentum fraction, while stands for the Bjorken scaling variable. 
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2 EXPERIMENTAL SET-UP 



Since the PGF process is dominated by quasi-real photoproduction (Q 2 — > 0), the perturbative scale for 
the selected events, p 2 , cannot be set to Q 2 as in the QCD analyses of inclusive data. Instead, this scale 
is chosen to be the transverse mass of the charmed quarks, p 2 = 4Mj = 4(m 2 + p\), where the D meson 
transverse momentum, pj, is defined with respect to the virtual photon. 




Figure 1: Photon-Gluon Fusion into a pair of charm quarks, cc. Symbols in parentheses denote four- 
vectors. 



The paper starts with a brief presentation of the experiment in Sec. [2] In Sec. [3] the data selection is 
reported in detail. The evaluation of the asymmetries and the corresponding results are described in Sec. 
[4] The determinations of the gluon polarisation (Ag/g) at leading (LO) and next-to-leading (NLO) QCD 
accuracies are presented in Sec. [5] Concluding remarks are given in Sec. [6] 

2 Experimental set-up 

The COMPASS spectrometer is a fixed target set-up situated at the M2 beam line of the CERN SPS 
using muon or hadron beams. For the present measurement, longitudinally polarised positive muons of 
160 GeV/c momentum were scattered off a large polarised solid state target. A detailed description of 
the set-up can be found in Ref. |[23l . 

The muons originate from the weak decay of 175 GeV/c pions and kaons produced by the 400 GeV/c 
SPS proton beam impinging on a primary beryllium target and are thus naturally polarised. The beam 
polarisation, P^, is about 0.8 at 160 GeV/c with a relative uncertainty of 5% Il24l . A beam intensity of 
about 4 • 10 7 muons/s was used, with the spill length between 4.8 s and 9.6 s for SPS cycles between 16.8 s 
and 48 s, respectively. The beam is focused onto the target centre with a spread of 7 mm (r.m.s.) and 
a momentum spread of 5% for the Gaussian core. The momentum of each incoming muon is measured 
with a precision better than 1% upstream of the experimental hall using a beam momentum station. 
Before the target, the trajectory of each beam particle is determined with an angular precision of 30 prad 
using a set of scintillating fibre and silicon detectors. 

The solid state target is housed in a large superconducting solenoid providing a field of 2.5 T with field 
uniformity, 8B/B, better than 10~ 4 . From 2002 to 2004, the angular acceptance was ±69 mrad at the 
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Table 1: Kinematic variables used to describe the muon-nucleon interaction, see Fig. 1 



Variable 


Symbol 


Definition 


Nucleon (muon) mass 


Mini?) 








Four vector of incoming (outgoing) muon 


k(k') 








Four vector of target nucleon 


V 

1 








Four vector of outgoing hadron final state 


PX 








Four vector of virtual photon 


q 


<? = 


k-kf 




Four vector of a final state hadron 


Ph 








Negative four momentum transfer squared 


Q 2 


Q 2 


= ~q 2 




(photon virtuality) 

Muon laboratory incident (final) energy 


E(E') 


E = 


(e> 

M V 


_p-k'\ 
M J 


Polar angle of kaon in the D° centre-of-mass 


d* 








relative to the D° laboratory momentum 










Virtual photon energy 


V 


v = 


V ■ a lnh 


-E' 


Bjorken scaling variable 


x B 


x B -- 


_ Q 2 

2pq 


Q 2 
2Mv 


Virtual photon fractional energy 


y 


y = 


p-q lab V 
p-k E 




Final state hadron fractional energy 


Zh 


Zh = 


P-Ph lab 

p-q 


Eh 
V 


Transverse momentum of D° meson 


D° 

Pt 








with respect to the virtual photon direction 










Energy of D° meson in laboratory 


E D o 









upstream edge and ±170 mrad at the downstream edge of the target material. From 2006 onwards, a 
new target magnet with a larger aperture solenoid was used [23 ] . It yields an angular acceptance of 
± 1 80 mrad for the upstream target edge resulting in a much improved hadron acceptance and matching 
the ±180 mrad acceptance of the spectrometer. 

The target material consisted of 6 LiD beads in 2002 to 2006 and NH3 beads in 2007, in a bath of 
3 He - He. The target was cooled down to a temperature below 100 mK by a 3 He - 4 He dilution re- 
frigerator. The target polarisation was accomplished using the method of dynamic nuclear polarisation 
(DNP) and measured continuously by a set of NMR coils surrounding the target material. The achieved 
polarisation, P t , was about 0.5 for deuterons ( 6 LiD) and 0.9 for protons (NH3) with a relative uncertainty 
of 5% and 2%, respectively. 

In 2002 to 2004, the target material was contained in two 60 cm long cells that were polarised in opposite 
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2 EXPERIMENTAL SET-UP 



directions. The polarisation was reversed 3 times per day by rotating the field of the target magnet. From 
2006 onwards, a three-cell target set-up was used with a central 60 cm long cell placed between two 
30 cm long ones. The material inside the central cell was polarised oppositely to that of the outer ones. 
The use of this new target arrangement allows for further reduction of the systematic uncertainty due to 
the variation of the spectrometer acceptance along the target, so that only one field rotation per day was 
performed. In order to minimise possible acceptance effects related to the orientation of the solenoid 
field, also the sign of the polarisation in each target cell was reversed several times per year by changing 
the DNP microwave frequencies. 

As not all nucleons in the target material are polarised, the so-called dilution factor, /, is introduced. 
It is expressed in terms of the numbers n A of nuclei with mass number A and the corresponding total 
(i.e. including radiative effects) spin-independent cross sections, a A \ per nucleon for all the elements 
involved: 

In the present analysis, the dilution factor is modified by a correction factor p = o^Jo^\ accounting for 
the dilution due to radiative events on unpolarised protons (deuterons) |[25l . A correction for polarisation 
of the deuteron in the 6 Li nucleus is also applied. 

The dilution factor depends on x#. At low x#, it is larger for events containing hadrons in the final state 
due to the absence of radiative elastic tails. Its values at medium xb for 6 LiD and NH3 are about 0.37 
and 0.14 with relative uncertainties of 2% and 1%, respectively . 

The two stages of the COMPASS set-up are open dipole spectrometers for large and small angle tracks, 
respectively. Each dipole magnet is surrounded by tracking detectors. COMPASS uses various types 
of them in order to match the expected particle flux at various locations in the spectrometer. In high- 
flux regions close to the beam, tracking is provided by arrays of scintillating fibers, silicon detectors, 
micromesh gaseous chambers and gas electron multiplier chambers. Further away from the beam, larger- 
area tracking devices as multiwire proportional chambers, drift chambers and straw detectors are used. 
In 2006 the tracking system in the first stage of the spectrometer was adapted to match the increased 
aperture of the superconducting solenoid. 

Muons are identified in large area tracking detectors and scintillators downstream of concrete or iron 
muon filters. Hadrons are detected by two scintillator-iron sandwich calorimeters installed in front of 
the muon filters. Electromagnetic lead glass calorimeters are placed in front of the hadron ones. The 
data recording system is activated by triggers indicating the presence of a scattered muon and/or energy 
deposited by hadrons in the calorimeters. Both inclusive and semi-inclusive triggers are used. In the 
former, the scattered muon is identified by coincident signals in the trigger hodoscopes, and in the latter 
the energy deposited in calorimeters is demanded in addition. Moreover, a calorimetric trigger with a 
high energy threshold is implemented to extend the acceptance. In order to suppress triggers due to halo 
muons, veto counters upstream of the target are used. The COMPASS trigger system covers a wide range 
in Q 2 , from quasi -real photoproduction to the deep inelastic region. 

For charged particle identification in the first stage of the set-up, a Ring Imaging Cherenkov detector 
(RICH) is installed ll26l . It is a gas RICH with a 3 m long C4F10 radiator. Two spherical mirror surfaces 
reflect and focus the Cherenkov photons on two sets of detectors situated above and below the acceptance 
of the tracking detectors, respectively. The photon detection uses MWPCs with segmented Csl photo- 
cathodes which detect photons in the UV region. In 2006, the central part of the RICH was upgraded 
replacing the MWPCs by multianode photomultiplier tubes yielding a considerably higher number of 
detected photons and a much faster response. For the outer parts, the readout electronics was refurbished 
allowing a significant reduction of the background. 

The particle identification procedure relies on a likelihood function based on information on the photons 
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detected in the RICH and associated with a charged particle trajectory. The likelihood function uses the 
photons of the signal and a theoretical expectation of their distribution, taking into account possible signal 
losses due to dead zones in the detector. For the description of the background photons, the experimental 
occupancy of the photon detectors is used. For each track, likelihood values are computed for different 
particle mass hypotheses and the background hypothesis. Identification of a pion (kaon) is possible for 
momenta between 2.5 GeV/c (9 GeV/c) and 50 GeV/c. 

The performance of the detectors as well as the stability of the reconstructed data was carefully monitored 
and all spills not fulfilling stability requirements were excluded from further analysis. Time intervals 
selected for asymmetry measurements correspond to periods of stable spectrometer performance. In 
total, data taking amounted to 48 weeks in the years 2002 to 2007. 

3 Data selection 

In order to extract information about the gluon polarisation, events with D mesons have to be selected 
from the data. This is accomplished by requiring every event containing an incoming and outgoing muon 
together with at least two outgoing charged tracks. Furthermore, only events with the incoming muon 
potentially crossing the whole target and with an interaction point (or 'vertex') within the target were 
retained. 

The direction of tracks reconstructed at an interaction point in the target is determined with a precision 
better than 0.2 mrad and the momentum resolution for charged tracks detected in the first (second) spec- 
trometer stage is about 1.2% (0.5%). The longitudinal vertex resolution varying from 5 mm to 25 mm 
along the target permits assigning each event to a particular target cell, i.e. to a specific target spin direc- 
tion. In Fig. [2] are shown the distributions of the reconstructed vertex position z vtx along the beam axis 
for events remaining after applying the aforementioned selection criteria. The relative increase of the 
number of events in the upstream and central target cells, seen in the right panel, reflects the acceptance 
increase due to the upgrade of the target magnet in 2006. 




Figure 2: Distribution of reconstructed vertex positions z vtx along the beam axis for the target with two 
(left) and three (right) cells. Dark horizontal bars at the bottom mark the target extensions, arrows denote 
the target polarisation directions. See text for details. 



Due to multiple Coulomb scattering in the solid state target, the spatial resolution of the vertex recon- 
struction is not sufficient to separate production and decay points of charmed mesons. As a result, such 
mesons can only be reconstructed using the invariant mass of their decay products. Their decay modes 
considered in this analysis are listed in Tablepl Those D° decays which involve the same set of final state 
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3 DATA SELECTION 



particles cannot be distinguished event by event. Therefore five independent data samples with different 
final states are defined for this analysis, see Table[3] In the four tagged samples, i.e. Djj^., Dk su1)7P D*^^, 
and D^ K7[Jl , the D° meson is assumed to originate from a D* decay into a D° meson and a slow pion, 

D* D°7T S . The kinematic selection criteria, which are tuned to reduce the combinatorial background 
without affecting the D° meson signal, are listed in Table |4] 

Table 2: Charmed D° meson decay modes, together with their branching ratios, considered in this anal- 
ysis. The charge conjugate (c.c.) final states from D° decays are also included. 



Reaction 
number 


D° decay mode 


1 


D° 


3.89% „_ j_ 

— > K 7T 


+ c.c. 


2 


D° 


13.9% T ,_ i n 


+ c.c. 


3 


D° 


8.09% T ,_ i i _ 


+ c.c. 



Table 3: Event samples used in the analysis. For each sample, the corresponding reactions from Table 
[2] are indicated. In the tagged samples, the D° is assumed to originate from D* decay and the final state 
of the D° decay is indicated by the subscript where 'K su t,' stands for a kaon with momentum below the 
RICH threshold. Throughout this paper, each sample will be referred to using the above notation. 





Untagged sample 


Tagged samples 


Sample 


D° 


u Ktc 


U K suh 7t 




n* 


Reaction number 


1 


1 


1 


2 


3 



Table 4: List of kinematic cuts used for each data sample. A D candidate is accepted if it fullfils 
all conditions in a corresponding column. Here AM =M^£ — —M % , where the superscript 'rec' 
denotes the reconstructed mass. 





Kinematic cut intervals 


Variables 


Untagged sample 


Tagged samples 


D° 




D* 

K<„/,Jr 


r>* 


(M--M D o)[MeV/c 2 ] 


[-400, +400] 


[-600, +600] 


[-400, +400] 


| cos 0*| 


<0.65 


<0.90 


<0.85 


Z D o 


[0.20, 0.85] 


[0.20, 0.85] 


[0.25, 0.85] 


[0.30, 0.85] 


PK [GeV/c] 


[9.5, 50] 


[9.5, 50] 


[2.5, 9.5] 


[9.5, 50] 


Pn [GeV/c] 


[7, 50] 


[2.5, 50] 


AM [MeV/c 2 ] 




[3.2,8.9] 


[4.0,7.5] 


P*[GeV/c] 




< 8 



Particles are identified using the RICH detector. Using the measured momentum of a charged particle 
and the distributions of Cherenkov photons, likelihood values for different mass hypotheses and for the 



9 



background hypothesis are computed. A particle is identified as kaon or pion if the likelihood value is 
larger than those for all remaining hypotheses. This procedure is very efficient in reducing the com- 
binatorial background of two particles other than % and K. A detailed description of the identification 
procedure is given in Ref. lf27l . 

The following selection criteria were applied to obtain the final event samples. The untagged sample 
contains events with Kn pairs in the reconstructed mass range given in Table [4j which do not stem 
from decays of reconstructed D* mesons, see Fig. [3] Due to large combinatorial background this sample 
requires more restrictive cuts for the identification of pion and kaon: a pion momentum above 7 GeV/c 
is required to avoid contamination from electrons. For the four tagged samples, a D* meson is selected 
by requiring the presence of a slow pion, p Ki < 8 GeV/c, in addition to a D° candidate. The presence 
of the slow pion permits the application of two additional cuts. The first one uses the RICH detector to 
reject electrons that mimic slow pion candidates and reduces the combinatorial background by a factor 
of two. The second one is a cut on the mass difference, AM = A/g^_ — — M n , where A/^L and 
A/JF£ are the reconstructed masses of the D* and the D° candidates, respectively. This mass difference 
can be measured with very good precision and thus the cut on AM results in a significant reduction of the 
combinatorial background in the tagged samples. 



xi<r 




100 200 300 400 

M Kit - M D „ ( MeV/c 2 ) 




M KlI - M D „ ( MeV/c^ ) 

Figure 3: Invariant mass spectra for the D^. and samples with the approximate number of D° 
mesons above background. 



In addition to the cuts described above, further kinematic cuts were applied to all samples. It is de- 
manded that |cos<9*| < 0.9 for the tagged Djj^. and D* « samples, |cos0*| < 0.65 for the sample 
and |cos0*| < 0.85 for the remaining samples. This cut suppresses mainly background events and im- 
proves the significance of the signal. Finally, all events have to satisfy a cut on z D o. Since a pair of 
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charmed quarks is produced in the centre-of-mass of the y*g system, each quark receives on average half 
of the virtual photon energy. Indeed, the measured z D o distribution and the one simulated assuming a 
pure PGF process (with parton showers included) are very similar and have a most probable value close 
to 0.5, see Fig. [4] This fact strongly supports the assumption on PGF dominance in charm production. 




Figure 4: Distribution of z D o for the Djj^ data sample (background is subtracted) and corresponding 
Monte Carlo events. D° mesons are selected in the ±80 MeV/c 2 mass window around the D° mass. 



Final mass spectra for the untagged and tagged samples, selected according to kinematic cuts listed in 
Table |4j are shown in Fig. [3] For the latter, a most pronounced signal is visible at the D° mass with a 
mass resolution of about 27 MeV/c 2 . In this mass spectrum, also a second structure is visible at about 
-250 MeV/c 2 which is due to events with D° — > K~n + n° decays with neutral pion not reconstructed 
in the analysis. Thus the Kn spectrum is shifted to lower mass as compared to D° — > K _ 7T + decays. 
The purity of this signal is much worse due to the non-reconstructed neutral pion. 

Further improvement of the significance of the signal is accomplished by applying the Neural Network 



method described in Sec. 4. 1.2 which leads to a considerable reduction of the combinatorial background 
in the tagged samples. The resulting mass spectrum for DJ^o is shown in Fig. |5j with an improvement 
of the signal strength by 15% while for the D° — > K~7i + the signal and the background are reduced 
in a similar way so that the significance of the signal stays unchanged. Therefore only the criteria from 
Table [4] are used to select the final Dj^ sample, see Fig. [3] 

Results on channels with a weaker D° signal like Djj^^ and % , are also shown in Fig. pi The sample 

limit of 5 



Dt _ contains events where the momentum of the kaon candidate is below the limit of 9 GeV/c for 
kaon identification by the RICH detector. Simulations using a Monte Carlo generator for heavy flavours, 
AROMA [28], and a full spectrometer description based on GEANT [29] have shown that about 30% of 
the kaons coming from D decays have their momenta below this RICH threshold. Therefore it is only 
required that those particles, K su b, are not identified as pions or electrons. 

In the case the two D° candidates are found in the same event, only one of them, chosen randomly, is 
considered in the analysis. If two channels contribute with a D° candidate to the same event, only one of 
them is accepted according to the following priority rule: Dj^^, Djj^ or D* D^, n (see Ref. 

E3). ^ ^ 

Distributions of xb, Q 2 and y variables for the D^ candidates from 2006, and from the D° signal region, 
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Figure 5: Invariant mass spectra for the Dj^o, ^>^ n7tn and D^ t)I samples. The purity of the samples 
was improved using the Neural Network. The approximate number of D° mesons above background is 
given. 



are presented in Fig. |6] For those events x# values range from about 10 5 to 0.1 with (jcb) = 0.004, Q 2 
values from 10~ 3 to 30 (GeV/c) 2 with (Q 2 ) = 0.6 (GeV/c) 2 , and y values from 0.1 to 1 with (y) = 0.63. 



4 Asymmetry evaluation 

In this section, we describe the determination of the virtual photon asymmetry for D° production, 
A rN = A rN ^ D x , from the event samples defined in Table ol The method is similar to the one used 
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Figure 6: A scatter plot of Q 2 vs xb and a distribution of y for the 2006 sample selected as in Fig. [4] 
except that the background is not subtracted. 



in our previous publication [19]. The asymmetry A^ N can be used in various ways to evaluate the gluon 
polarisation (Ag/g) at LO or NLO QCD accuracy. 

4.1 Analysis method 
4.1.1 Asymmetries 

The number of events collected in a given target cell and time interval is given by: 
d k N 



■ a$n(s + b) 
dmdX 



i + w ( 4^ n ^' d ° x + ±^ B ) 



(1) 



Here, ^ N ^^' D ° X is the longitudinal double spin asymmetry of the differential cross section for events 
with a D° or D° in the final state, and Ab is the corresponding asymmetry originating from background 
events. Furthermore, m = M^ n (or m = M^nn) and the symbol X denotes a set of k — 1 kinematic 
variables describing an event (pj , £" D o, Q 2 , y, ...), while a, <p and n are the spectrometer acceptance, the 
incident muon flux integrated over the time interval, and the number of target nucleons respectively. The 
differential unpolarised cross sections for signal and background, folded with the experimental resolution 
as a function of m and X, are represented by s = s(m,X) and b = b(m,X) respectively. The ratio s/(s + 
b) will be called 'signal purity' and the ratio b/(s + b) 'background purity'. The information on the 
gluon polarisation is contained in the virtual photon asymmetry A rN = A^ N ^^' D X /D. Similarly, the 
background asymmetry can be written as Ag N = A^/D. Here, D is the so-called depolarisation factor 
accounting for the polarisation transfer from lepton to virtual photon: 



D 



y 


2-y 






Q 2 






Q 2 


)+2(l 


-y) 



(2) 



A straightforward way to extract A rN would be the following: Eq. {!]) is integrated over the variables X 
to obtain the number of events in both spin configurations as a function of the invariant mass m. Next, 
the event number asymmetry in the D signal region is extracted, and a possible background asymmetry, 
determined from the asymmetries in sidebands to the left and right from the signal region, is subtracted. 

For this analysis, however, we choose the method of event weighting, which is advantageous in terms 
of statistical precision. Compared to previous COMPASS analyses where weighting procedures were 
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applied (3), here the weighting procedure is extended to determine the background asymmetry Ag N si- 
multaneously with A rN [30]. In order to achieve this, every event is weighted once with a signal weight, 
ws, and once with a background weight, wb- 

ws = PufD^r, (3) 
s + b 

wb = Pnf D -^T- (4) 

s + b 

Except for the target polarisation P t , these weights are the prefactors of the asymmetries A rN and Ag N , see 
Eq. ([T]). The target polarisation is not included in the weights because its time dependence would lead 
to an increase in false asymmetries. The signal and background purities are included in the respective 
weights. This procedure leads to the highest possible statistical precision which would also be obtained 
in the unbinned maximum likelihood method ll30l |3T1 . Note that the unbinned maximum likelihood 
method, cannot be applied here because the acceptance and flux factors in Eq. ([T]) are not known with 
sufficient precision, only their ratios for different spin states and target cells are known. These ratios will 
be used for the extraction of A yN and Ag . 

The asymmetry A rN is extracted in bins of D° transverse momentum with respect to the virtual photon, 
Pj , and D° energy in the laboratory system, E D o. These variables and their binning were chosen to 
minimise the influence of the experimental acceptance on the asymmetry. As A rN does not contain the 
depolarisation factor D, its remaining dependence on the inclusive variables y and Q 2 is very weak. The 
expectation value of the sum of signal weights is obtained as: 




w s (X,m) 



d k N t 
dmdX' 



dmdX 



as,, 



The symbols used are defined as: 



as.f = J ws a t §t n t (s + b) dm dX , 
/ rjwa t (j)tn t (s + b)dmdX 



/ wa t (j) t n t (s + b)dmdX 



(5) 

(6) 
(V) 



with j8 s = w s P t , j8 B = w B Pt, tj e jSs^B^VJT and 

w € [ws,wb,wsj8s,wsj3b]. The index t denotes 
the target cell before (t = u,d) or after (t = u',d') spin rotation and N t is the number of events observed 
in a given target celQ An equation analogous to Eq. §5§ holds for the sum of background weights 
(£i=i w bj) » with analogous definition of symbols. In total, eight equations similar to Eq. (Dj) are obtained 
for every (pj° ,E D o) bin: for the signal and background weights in two target cells and for two spin 
configurations. These eight equations contain 12 unknowns: (A rN ) H , sjSs , (Ag N ) H , sftj , (A rN ) M , B& , (A B fi ) WB ^ B , 
four acceptance factors as,t and four acceptance factors otB.t- 

The expectation values of the sum of weights on the left hand side of Eq. (|5]> are identified with the 
measured sums of weights. In order to extract A rN and Ag N from the measured sums of weights one 
proceeds as follows. The factors (J3s,b)w s ,w b 316 evaluated from data, e.g. 



w s 



(8) 



2 In 2006-2007 d and u stand for the central target cell and the sum of the outer ones, respectively. 
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The expectation values appearing in Eq. (|7]) contain the spin-averaged cross section while the sum over 
events in Eq. ([8]> used to evaluate these expectation values runs over the spin-dependent events. This has 
a negligible effect on the result because the raw asymmetry P t P^fD(A rN } is very small. This smallness 
makes sure that neither the result nor its statistical error are sensitive to the fact that the same data which 
are used to determine the asymmetries are also used to evaluate the expectation values above. 

The acceptance factors as,t and aB,t cannot be determined with sufficient precision to extract A rN and 
A^ 1 directly from the set of eight equations. By assuming that for both signal and background possible 
acceptance variations affect the upstream and downstream cells in the same way, i.e. a u /a d = a u /a d , 
the number of unknowns is reduced to ten. With an extra, much weaker assumption that signal and 
background events from the same target cell are affected in the same way by the acceptance variations, 
one arrives at a system of eight equations with nine unknowns. Possible deviations from the above 
assumptions may generate false asymmetries which are included in the systematic uncertainty, see Sec. 
EQ 

The number of unknowns is reduced to seven with two additional assumptions: 

(A? N ) R =(A^) H =A? N and «\ = « \ = Af, (9) 

which are satisfied for constant values of A rN and Al N in a given bin. The uncertainty on the gluon 



polarisation introduced by this assumption will be discussed in Sec. 5.2 Using the set of eight equations, 
the asymmetry A rN and the background asymmetry Ag N are determined simultaneously by a standard 
least square minimisation procedure, which takes into account the statistical correlation between £ws 
and £wb in the same target cell. The correlation factor cov(£ws,Lwb) is given in Ref. [31]. The 
analysis is performed independently for each (pj ,E D a) bin and D meson decay channel. 



For determinations of average values of kinematic variables in each (pj ,E D o) bin, a weight equal to w| 
is used, in accordance with Eq. ([5]). 

4.1.2 The signal purity 

The signal purity s/{s + b) can be extracted from a fit to the invariant mass distribution of D° candidates. 
It depends on kinematic variables, for instance it is large at high transverse momenta p^ of the Kit 
system and small at low pj . In order to implement the kinematic dependence of the signal purity in the 
weights given by Eqs (|3]|4]>, one would naively proceed by performing fits to the corresponding invariant 
mass distributions in bins of kinematic variables. This procedure is not feasible in our case because of 
limited statistics. Instead, in this analysis a classification based on a Neural Network is employed E71 . 

Here, the aim of the Neural Network is to distinguish signal from background events using only data. 
The network consists of information processors (neurons), which are interconnected and organized into 
layers. The external information fed into the input layer is processed in the hidden layers and the result 
produced by the output layer is a classification of the event by the network. In the present case, the 
input layer contains a set of kinematic observables: ratios of RICH likelihoods, cos0*, z.d and kaon 
momentum. There are two hidden layers and the number of neurons in them varies during the training 
process (dynamic network). For each event, the network tunes the strength of each variable-neuron and 
neuron-neuron connection. The strengths are obtained by minimising the squared deviation between the 
expected output and the actual Neural Network prediction. This training process stops when the deviation 
reaches a stable minimum ll32l . 

For each event sample (see Table[3]), two data sets are used as inputs to the Neural Network. The first one 
contains the D° signal and the combinatorial background events. These events are called 'good charge 
combination' ones (gcc) referring to the charges of particles from D° decays, and they are selected as 
described in Sec. [3] The second set, the 'wrong charge combination' events (wcc), is selected in a 
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M Kn - M D „ ( MeV/c 2 ) 

Figure 7: Example of the Kit invariant mass spectrum for the Djj^ sample with good (gcc) and wrong 
(wcc) combination of pion and kaon charge signs. See text for details. Data were collected in 2007 with 
the proton (NH3) target. 

similar way except that the sum of charges of corresponding particles should not be zero. It contains 
only background events and is used as a background model (see Fig. [7]). The Neural Network performs 
a multi-dimensional comparison of gcc and wcc events in a ±40 MeV/c 2 mass window around the D° 
mas^] Within the gcc set, signal events are distinguished from combinatorial background by exploiting 
differences between the gcc and wcc sets in the shapes of distributions of kinematic variables as well as 
multi-dimensional correlations between them. An example of a properly chosen variable for the network 
is cos 8*, as shown in Fig. [8] The reconstructed mass cannot be used because it would enhance the 
probability of a background event in the signal region to be a true D meson. 

The network classifies all the gcc events according to their similarity in kinematics with respect to the 
wcc ones, and to each event it assigns a probability of being a signal. A probability of 0.5 is assigned to 
indistinguishable events. If the network is trained with proper input samples, i.e. a correct background 
model and a sufficiently strong signal, the network output, [s/(.s + &)] NN , can be directly interpreted as 
an estimate of the signal purity in the corresponding mass window. This is the so called 'pure' Neural 
Network method, applicable to the Dt „» and D^ nnn samples collected in 2004-2007, where event 
statistics and signal purities are large. 

The mass dependence of signal and background strengths, s(m) and b(m), which cannot be obtained 
from the Neural Network in an unbiased way, is determined from a fit to the mass spectra in bins of 
[s/(s + &)] NN . In order to describe the signal a Gaussian distribution is used for all samples, while for 
the background the following fitting functions are employed: two exponential distributions for the 
channel and one exponential for the D* tagged channels. An exception is the sample, for which a 

second degree polynomial is used. From those fits, corrections X to the signal purity are obtained in the 
mass windows defined above: 



A mass window of ±30 MeV/c 2 is used for the sample and ±40 MeV/c 2 around -250 MeV/c 2 for the sample D* 
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Figure 8: Example of the distribution of | cos 6*\ (polar angle of kaon in the D° centre-of-mass relative 
to the D° momentum) in the D° meson rest frame for the gcc and wcc events (Djj^. sample, 2006 data). 
Top: region of the D° signal, bottom: outside the D° signal. 



so that 



, , b(m)dm 

/ s(m)dm 



Xs(m) 



s + b Xs(m) +b(m) 



(10) 



(11) 



The fit of the invariant mass spectra in bins of the NN signal purity can also be used to validate the 
classification obtained by the Neural Network. For each bin, the signal purity is determined from an 
integration of the signal and background fits over the mass windows used. Good agreement between 
signal purities from the NN and the fit is found for all samples, which confirms that the Neural Network 
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does not introduce any bias in the analysis. As an illustration, the mass spectra in bins of the NN signal 
purity together with a comparison of the two signal purities are shown for the Dj^ sample in Fig. [9] 
The signal purity clearly increases with increasing [s/ (s + &)]nn- Equally good agreement is found when 
comparing s/(s + b) between NN and data in bins of (pj , E D o). 
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Figure 9: The Kti invariant mass spectra in bins of the NN signal purity [s/ (s+&)]nn for the Djj^. sample. 
The last panel shows a comparison of the two purities, [s/ {s + b)]^ and [s/ (s + b)]^ (see text for details). 
Curves show the background component of the invariant mass tits described in the text. The significance 
of the signal is shown as ratio s/b. 



The signal purity can be parameterised in various ways provided it correctly reproduces the data. Several 
parameterisations were found to indeed yield asymmetries compatible within statistical uncertainties. 
In order to achieve the statistically most precise result on the gluon polarisation, we chose the method 
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described in Sec. 5.2 As a special feature in this method the signal weight contains the product of the 
signal purity and all, which is the partonic asymmetry for the PGF process. Observation shows that 
£?ll is strongly anticorrelated with the signal purity. Therefore a parameterisation of ajx, validated as 



described in Ref. [27], was additionally included in the training of the NN. 

For the low purity sample collected in 2002-2007, and for all samples collected in 2002 and 2003, 
the extraction of the signal purities from the network is more complicated since the anticorrelation men- 
tioned above cannot be accounted for because of weak signals. Therefore a 'hybrid' method is employed. 
Similar to the method used in Refs | :T9l[33l . this approach uses fits to the mass spectra which are sam- 
pled in bins of two variables, NN signal purity (from a parameterisation without ajx in the training) and 
fP^LL- The former sorts the events according to their similar kinematic dependences, while the latter is 
used to ensure the anticorrelation between «ll and the signal purity. The signal and background distribu- 
tions belonging to the mass spectra sampled in the bins are fitted by the same fitting functions as defined 
above for describing the mass dependence of signal and background. Integrating the fits within the same 
mass windows as used for the NN training procedure yields the signal purities extracted from the fit. For 
each of the two variables, a function is built using linear interpolations between the fit results. An itera- 
tive procedure is used to obtain a stable result on these two functions simultaneously, and thereafter the 
correction X is applied to the signal purity. Due to the statistical limitations, only one parameterisation 
was built for each decay channel and year. 

As the hybrid method can be used for both, low purity and high purity channels, it was decided to use it 
for all parametrisations of signal purities. Although more complex than the pure NN method, the hybrid 
method results in a comparable statistical precision. 



4.2 Results on asymmetries 

The asymmetry A rN is extracted simultaneously with the background asymmetry Ag N for each bin, chan- 
nel, and year of data taking except for low purity channels where some data taking years of the same 
target set-up are merged. Final results sorted by D° decay mode are shown in Tables [5]— |7j where A^ 
is given in each (pj ,E D o) bin together with average values of kinematic variables. All averages are 
calculated with the weight w|. The muon-nucleon asymmetry A^ N ^' D x can be obtained from A yN by 
multiplying it by (D), which is a function of (y) and (Q 2 ) and is also given in Tables|5]-|7] 

As the extraction of A rN is performed based on event weights, uncertainties introduced into the determi- 
nation of the asymmetry by each contribution to the weight have to be accounted for. The uncertainty 
of A rN is acquired from a spread of weighting factors w, Eq. which is obtained by comparing the 
default analysis with weight wo to other analyses with different weight, w. The expectation value of the 
weighted asymmetry is (A^) = (tvtvo) / (wq)A 7N , as shown in Ref. 11271 . The spread of (wwq) / (wq) gives 
the relative systematic uncertainty of A rN . 

The major sources of systematic uncertainties in the measurement of A rN are discussed below. The con- 
tributions from Pp, P x and / are taken conservatively as 5%, 5% and 2% respectively, for both deuteron 
and proton targets. In order to study the contribution of s/{s + b) to the systematic uncertainty, three 
tests are performed. In one of them, different fitting functions are used for the functional form of the 
background. In the other two, different mass windows are investigated concerning the Neural Network 
parameterisation and the choice of the binning used in the reconstruction of the D° spectra. Note that 
each of these tests leads to several new values of s/(s + b) and, consequently, new weights w are ob- 
tained. The resulting spread of weights is computed for each year of data taking, each sample and each 
bin with respect to the default weight wq. Thereafter, the weighted average of all spreads is determined 
separately for each of the three systematic tests considered. The combined uncertainty on s/ (s + b) is 
obtained from a quadratic superposition of these three uncertainties. The resulting average value over all 
bins is 7% of the measured asymmetries. 



4.2 Results on asymmetries 
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The contribution of D to the uncertainty of A rN is obtained as follows. According to the experimental 
uncertainty of 1% in the measurement of the mean value of the scattered muon momentum, shifted values 
of y are calculated for every event. Thereafter, new values of D are computed from Eq. Q. The resulting 
spread of (wwo) / (wq) gives a systematic uncertainty of 1.6%. 

Systematic uncertainties of A rN arising from false asymmetries and from the assumptions specified in 
Eq. Q can be best estimated using the statistically optimised method (see Sec. 5.2 1. First, they are 
determined for the gluon polarisation Ag/g, and then they are translated to A rN in bins of (pj ,E D o) 
employing (a^/D). Averaging over all bins, the resulting absolute values of the uncertainties due to 
false asymmetries and the assumptions in Eq. Q are 0.022 and 0.007, respectively. 

The total systematic uncertainties of A 7N , as given in Tables|5]-|7] are obtained by adding all contributions 
in quadrature. 
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5 DETERMINATION OF THE GL UON POLARISATION 



5 Determination of the gluon polarisation 

In this Section we present the results of our measurement of the gluon polarisation. The extraction of 



Ag/g from A^ N at LO QCD accuracy is discussed in Sec. 5.1 The LO determination of the gluon 



polarisation by a statistically optimised method is described in Sec. 5.2 The extraction of Ag/g from 



A rN at NLO accuracy is presented in Sec. 5.3 



This analysis neglects any contribution from 'intrinsic charm', i.e. nonperturbative charm quark or 
charmed hadron components of the nucleon wave function. Such contributions, estimated to be < 1% 
04ll35l . are fundamentally different from the perturbative splitting of a gluon into a cc pair; the latter 
decreases strongly with xb- In the EMC measurement of the charm component in the nucleon structure 
function F 2 CC ll20l . a possible intrinsic charm contribution of about 1% at xb ~ 0.4 could not be excluded 
||20l l2TTl . Up to now, the estimates of Refs ||34l 1351 cannot be experimentally verified due to the poor 
statistics of the EMC measurement at large xb, too low values of xb in the HERA F 2 CC measurements 
|[36l . and kinematic acceptance limited to the region x% < 0.1 for open charm production in COMPASS. 



The contribution of resolved-photon interactions was estimated using the RAPGAP generator 11371 and 
found to be negligible in our kinematic domain. 

5.1 Leading Order results from the asymmetries 

The information on the gluon polarisation contained in A^ N can be decomposed at LO accuracy as 

= DA^ = aLL — , (12) 
g 

assuming photon-gluon fusion as the underlying partonic process. Here «ll is the analysing power (also 
called 'partonic asymmetry') of the }ig — > fi'cc process. 

The analysing power £?ll depends on partonic kinematics. It is not accessible experimentally on an 
event-by-event basis. It is obtained using the Monte Carlo generator AROMA [28] in leading order 
QCD approximation, i.e. with parton showers switched off. The generated D events are processed with 
GEANT |[29l to simulate the full response of the COMPASS spectrometer, and then are reconstructed 
with the same analysis chain as used for real events. In order to provide £?ll values for real data, a Neural 
Network with the same architecture as described in Sec. [4] is used to parameterise the generated «ll 
in terms of measured kinematic variables X. Here the input layer contains the following observables: 
X = {Q 2 , y, x, pj , E D o}. As a result, a^(X) is obtained for real data on an event-by-event basis. 

Contrary to the parameterisation of s/(s + b), the Neural Network predicts values for c?ll based on 
event kinematics. For each generated event, the network tunes the strength of each variable-neuron and 
neuron-neuron connection. The strengths are obtained by minimising the squared deviation between the 
expected output, i.e. the generated «ll> and the actual Neural Network prediction for «ll based on X. 
This training process stops when the deviation between generated and parameterised cjll reaches a stable 
minimum. Six separate au^{X) parameterisations were built: for three D meson decay channels (K7T, 
K7r7r°, Ktititi), each for two experimental configurations (2002-2004 and 2006-2007). The correlation 
achieved between the generated and the parametrised analysing powers is 77% for the Dj^o channel 
and 82% for the remaining channels. The trained network is applied to real data. 



Knowing that A^ N = DA rN , the determination of the gluon polarisation from Eq. ( 12 ) is straightforward. 
The extraction of Ag/g from A rN in bins of (pj , £" D o) is performed using the values of (aix/D) shown 
in Table [8] 

The gluon polarisation in LO QCD, obtained from A rN amounts to 

'As x 

1 0. 10 ± 0.22 (stat.)± 0.09 (syst.) (13) 



5. 1 Leading Order results from the asymmetries 
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in the range of 0.06 < x < 0.22 with a weighted (x) « 0.11, and a scale (/I 2 ) « 13 (GeV/c) 2 . The range 
of x is determined by the r.m.s. value of a Gaussian distribution in log 10 x. Assuming that Ag/g(x) is 
approximately a linear function of x in the range covered by the present data, the above result corresponds 
to the gluon polarisation Ag/g at the value (jc). 

The sources of systematic uncertainties in the measurement of Ag/g are listed in Table [9j The contri 



butions from P^, P t , f, s/{s + b) and D are the same as discussed in Sec. 4.2 Contributions from false 



asymmetries and from the assumption given in Eq. ([9]), the same as in the statistically optimised method, 



are discussed in Sec. 5.2 In order to estimate the influence of the simulation parameters on the determi- 
nation of <3ll> Monte Carlo samples with different parameter sets are generated and the analysing power 
is recalculated. In these parameter sets, the mass of the charm quark is varied between 1.3 GeV/c 2 and 
1.6 GeV/c 2 , and the parton distribution functions as well as the factorisation scale is varied by a factor 
of eight. From each of these systematic tests a new value of {ai^/D) is obtained, and thereafter Ag/g is 
recalculated for each (pj , E D o) bin by dividing A rN by (au^/D). The systematic uncertainty in each bin 
is determined from the average spread of Ag/g compared to the result of the default analysis. The value 
for the systematic uncertainty of gluon polarisation is obtained as a weighted average of the systematic 
uncertainty in each bin. The relative uncertainty introduced by cjll alone is 15%. 



The final systematic uncertainty of {Ag/g) is obtained as a quadratic sum of all contiibutions. 
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Table 8: The average LO photon-gluon asymmetries, (au^/D), in bins of (pj , E D a) for each D° decay 
mode studied in the analysis. The averages use a\^/D from data events, obtained from the Neural 
Network parameterisation; they are weighted with w\. 



Bin limits 


Photon-gluon asymmetries 


D° 

Pr 


£" D o 


DL- d *kk and D* KsubJl samples 


D K**° sam P le 


d kxkt: sam P!e 


(GeV/c) 


(GeV) 


combined 






0-0.3 


0-30 


0.65 


0.62 


0.64 


0-0.3 


30-50 


0.68 


0.65 


0.63 


0-0.3 


> 50 


0.76 


0.74 


0.74 


0.3-0.7 


0-30 


0.46 


0.42 


0.38 


0.3-0.7 


30-50 


0.50 


0.46 


0.41 


0.3-0.7 


> 50 


0.56 


0.53 


0.52 


7 1 




96 


1Q 




0.7-1 


30-50 


0.26 


0.21 


0.25 


0.7-1 


>50 


0.29 


0.26 


0.30 


1-1.5 


0-30 


0.00 


-0.06 


0.02 


1-1.5 


30-50 


0.01 


-0.05 


0.04 


1-1.5 


>50 


0.05 


-0.02 


0.08 


> 1.5 


0-30 


-0.23 


-0.29 


-0.26 


> 1.5 


30-50 


-0.26 


-0.31 


-0.23 


> 1.5 


>50 


-0.27 


-0.31 


-0.22 



Table 9: Contributions to the systematic uncertainty of (Ag/g) LO obtained from A rN . 



Source 




Source 




Beam polarisation 


0.005 


s/(s + b) 


0.007 


Target polarisation P t 


0.005 


False asymmetry 


0.080 


Dilution factor / 


0.002 




0.015 


Assumption, Eq. (|9j) 


0.025 


Depolarisation factor D 


0.002 


Total uncertainty 0.086 
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5 DETERMINATION OF THE GL UON POLARISATION 



5.2 Statistically optimised determination of the gluon polarisation at LO 



The data described in Sec. [5] allow for the determination of (Ag/g) in a different, statistically optimised 
way. Practically it means that the gluon polarisation is obtained by replacing the factor D with a^L in 
the definition of ws- ws = PufctuJ / (s ± b). The use of this weight allows us to reproduce the results 
on (Ag/g) obtained from A yN ((/7°°), (E D o)) with about 6% gain in the statistical precision. This gain is 
due to a wide range of «ll values but the observed (anti)correlation between the signal purity and the 
parameterised a^L has to be accounted for in the parameterisation of s/(s + b). This fact is crucial to 
obtain an unbiased result of (Ag/g) with this statistically optimised method. 

Values for (Ag/g) and the background asymmetry (Ag N ) were obtained for each of the 48 weeks of data 



taking and separately for each of the five event samples. The results shown in Table 10 are the weighted 
means of those values. 



Table 10: Results for (Ag/g) and (Ag N ) for each data sample. Errors are statistical. 





u Ktc 




u Kra 


D* 


D° 


(Ag/g) 


-0.192 ±0.305 


-0.414±0.575 


0.614 ±0.667 


0.497 ±0.995 


0.020 ±0.415 


O 


+0.019 ±0.029 


±0.051 ±0.035 


±0.004 ±0.036 


±0.004 ±0.047 


-0.005 ±0.004 



The value of the gluon polarisation is obtained as the weighted mean of the five results shown in Table 
[10] and amounts to 

A " ' y -0.06 ±0.21 (stat.) ±0.08 (syst.) (14) 



in the range of 0.06 < x < 0.22 with a weighted (x) w 0.11, and a scale (/I 2 ) « 13 (GeV/c) 



The major contributions for the systematic uncertainty given in Eq. ( [T4] ) are presented in Table 1 1 
They were estimated as follows. In addition to P^, P t , f and s/(s + b), the uncertainty of «ll was also 
determined from the spread of weights (wwq) / (\Vq) (where wq stands for the default analysis). The use 



of different sets of parameters as described in Sec. 5.1 gives rise to a relative systematic uncertainty of 



9% of the (Ag/g) value originating from <3ll- The relative systematic uncertainty introduced by s/ (s + b) 
is 7%. 

In order to study the influence of false asymmetries, the Djj^ sample was divided into two sub-samples 
using criteria related to the experimental apparatus, e.g. the slow pion going to the left or to the right side 
of the incoming muon. The resulting asymmetries were found to be compatible within their statistical 
accuracies, i.e. no false asymmetries were observed. An upper limit of the contribution of time dependent 
acceptance effects to the systematic uncertainty was derived from the dispersion of (Ag/g) and (Ag N ) in 
the 48 weeks of data taking. The study was performed using the background asymmetry, profiting from 
the large statistics. Then the obtained results were translated to (Ag/g) using the method described in 
Ref. (27). An uncertainty of 0.024 was obtained assuming that possible detector instabilities are similar 
for background and signal events. Notice that the same assumption was used to reduce the number of 
unknowns in Eq. Q from 12 to 9. Therefore, a more conservative approach is taken: the double ratio of 
acceptances for the signal, a$ • a| / a| • , is assumed to be uncorrected with the corresponding one 
for the background events. The combination of these two cases leads to an upper limit of 0.08 for the 
possible contribution of false asymmetries. This contribution is also used in Sections |4~2~| and [BTTj 

An uncertainty originating from the assumption for Ag/g, analogous to that given in Eq. ^ and valid for 
Ag/g constant in the measured interval, 0.06 < x <0.22, is estimated as follows. First a pair of extremal 
values of Ag/g in that interval is selected as those given by the COMPASS NLO QCD fit with Ag > 0, 
see Sec. |5.3| This fit is chosen to maximise a potential influence of the above assumption. Next, a 
difference between these two Ag/g values is used as a bias in the system of equations from which Ag/g 
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Table 11: Contributions to the systematic uncertainty of (Ag/g) LO obtained in a statistically optimised 
method. 



Source 




Source 




Beam polarisation 


0.003 


s/(s + b) 


0.004 


Target polarisation P t 


0.003 




0.005 


Dilution factor / 


0.001 


False asymmetry 


0.080 


Assumption, Eq. 


0.025 






Total uncertainty 0.084 



and Ag N are determined. The bias is added to all (Ag/g)^^ terms while the (Ag/g) ws p s ones are left 
unchanged. The Ag/g interval resulting from the equations gives a relative systematic uncertainty due 
to the assumption. Possible variations of Ag N are studied in a similar way using a parameterisation of 
the inclusive asymmetry, Af, lf38l . The systematic uncertainty on Ag/g is taken as the largest difference 
between the result obtained in the default analysis and results of those tests. 



The result given in Eq. ( \hty is shown in Fig. [TO] together with a compilation of other LO gluon polari- 
sation measurements from high-/? T hadron production by COMPASS 03] [HI, SMC [H3| and HERMES 
lfT4l . The present measurement is at a scale of about 13 (GeV/c) 2 while other measurements are at 3 
(GeV/c) 2 . 



O) - 8r 
< 

0.4 - 
0.2 - 



■ COMPASS, high p t , Q 2 >1 (GeV/c) 2 , 2002-2006 
* COMPASS, high p , Q 2 <1 (GeV/c) 2 , 2002-2003 

COMPASS, open charm, 2002-2007 
□ SMC, high p t , Q 2 >1 (GeV/c) 2 

HERMES, high p , all Q 2 



o ■ 

-0.2 ■ 
-0.4 ■ 
-0.6 ■ 



-0.8 



10" : 



J I I 



10" 



Figure 10: A compilation of gluon polarisation measurements from open charm and high-^j hadron 
production. The star denotes a result of the present, open charm analysis, Eq. (14), obtained at LO 
accuracy, 2002-2007 data and all values of Q 2 . Full squares denote a COMPASS result JTBl for high-px 
hadron production on 2002-2006 data, for Q 2 > 1 (GeV/c) 2 while a full circle corresponds to 2002- 
2003 data and Q 2 < 1 (GeV/c) 2 [15]. The empty square shows the SMC measurement 0J] for Q 2 > 1 
(GeV/c) 2 and the empty triangle the HERMES result [14] obtained for all values of Q 2 . The horizontal 
bars mark the range in x for each measurement, the vertical ones give the statistical precision or the total 
uncertainties. 
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5.3 Next to Leading Order results 



The extraction of the gluon polarisation as described in Sees 5.1 - 5.2 was performed at LO where the 
only process leading to open-charm production is PGF. This requires knowledge of the analysing power 
oll and the signal strength on a bin-by-bin or event-by-event basis. Only combinatorial background was 
considered in the LO analysis. 

In this section a brief outline of a method of computing the NLO QCD corrections to the £?ll calculation 
in our analysis is given, see also Ref. |[39l . Examples of the NLO processes contributing to the muo- 
production of the cc pair are shown in FigjTT] Apart from the NLO corrections to the PGF mechanism, 
FigrjTJa-c, there exists yet other NLO contributions to muoproduction of open-charm, initiated by light 
quarks; as an example a process where a gluon emitted by a light quark creates the cc pair is shown 
in Fig. [TT] l Such processes do not probe the gluons inside the nucleon albeit they contribute to the D 
meson signal. Therefore, in the extraction of the gluon polarisation at NLO accuracy from the signal 
asymmetries, a correction term has to be taken into account: 



A 7N 



all Ag 
D 8 



Here, all is calculated at NLO accuracy and is different from the corresponding one at LO. 





(a) 



(b) 





(c) 



(d) 



Figure 1 1 : Examples of the NLO processes contributing to the muoproduction of the cc pair: a) virtual 
correction, b), c) gluon bremsstrahlung, d) light quark background 



The QCD calculations at NLO accuracy for spin averaged [40] as well as polarisation dependent cross 
sections for open charm production [41] are available only at the photoproduction limit, i.e. for Q 2 = 0. 
They are used in our analysis to estimate the value of the NLO corrections to c?ll and the light quark 
contribution. The average value of Q 2 in the kinematic region of our measurement is about 0.6 (GeV/c) 2 . 
It was confirmed by a direct check at LO accuracy that the Q 2 = limit used in the calculation is a very 
good approximation in our kinematic domain. 

Note that only one D meson is registered in the COMPASS data, while the second charm particle is 
unobserved. Also the NLO calculations of Refs BUI l4lTl represent integrated cross sections for a single 
charm quark (meson) observed in the final state. The cross sections are integrated over the kinematic 
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variables of the 'unobserved' second charm quark and radiated hard gluon, present in NLO processes. 
The limits of the integration depend on the available phase-space, left for 'unobserved' partons which is 
determined by x and the kinematics of the D-tagged c quark. 

In order to obtain <3ll in NLO accuracy on an event-by-event basis, the AROMA generator (which is 
based on a LO matrix element) with parton showers included is used, followed by a full simulation of the 
detector. In this way an approximation to the phase space needed for NLO QCD corrections is provided. 

For every simulated event, the upper limit of the integration over the energy of the unobserved gluon 
in the NLO emission process, y*g — > ccg, is obtained from the partonic Mandelstam variables, s and 
t. Both variables are calculated from the kinematics of simulated events. In particular s is determined 
from x, x# and Q 2 while t is related to the kinematics of the D-tagged charm quark. The integration 
over unobserved NLO real gluon emission reduces a differential cross section for a three-body final state 
(ccg) to that for a two-body one (cc), which has to be added to the LO cross section (cc, PGF) and the 



two-body virtual corrections (see e.g. diagram (a) in Fig. 11 ). In this way a correct infra-red divergence 
cancellation is achieved BP . The semi-inclusive partonic cross section at NLO accuracy is calculated 
on an event-by-event basis for both spin averaged and spin dependent case using formulae of Ref. PT1 . 
and consequently a^L at NLO accuracy is obtained. The same procedure is applied for the correction 
originating from a light quark. It should be stressed that in this method of <3ll estimation at NLO, only 
the values of s and t are taken from AROMA simulated events. 

To obtain the gluon polarisation at NLO accuracy, the measured asymmetry for D meson production has 
to be combined with the aix calculated at NLO. The kinematic variables pj , £" D o and gluon momen- 
tum fraction x define the total energy of all particles produced in the final state of the partonic process, 
including unobserved gluons emitted at NLO accuracy, y*g — > ccg. However, in o^l calculations, sim- 
ulated events at given x values are used and the integration over energy of the unobserved gluon is 
performed to obtain divergence-free a^L, which depends on § and t only. Therefore to be consistent with 
the NLO method of calculating oll, the asymmetry measured in (pj , E D o) intervals is binned into five 
one-dimensional intervals of pj° only. 



Table 12: The average values of the (o"ll/-D) and A corr , at NLO and in bins of pj° for each D° decay 
mode studied in the analysis. 



D° 

Pt 


n° d* 

U K7l> U Kk 


and b7t samples 


D k™° sam P le 


D M Sam P le 


(GeV/c) 




combined 












{cill/D) 


^corr 


{oll/D) 


Acorr 


(ull/D) 


Acorr 


0.0-0.3 


-0.130 


0.001 


-0.127 


0.002 


-0.097 


0.000 


0.3-0.7 


-0.241 


0.003 


-0.263 


0.003 


-0.240 


0.001 


0.7-1.0 


-0.419 


0.005 


-0.460 


0.004 


-0.404 


0.002 


1.0-1.5 


-0.574 


0.008 


-0.607 


0.008 


-0.572 


0.006 


> 1.5 


-0.679 


0.027 


-0.710 


0.020 


-0.719 


0.021 



In each p® bin, the weighted averages of a^/D and A CO rr are calculated, Table 12 and the gluon 



polarisation is evaluated from the A rN ^ D ° x asymmetries. The NLO light quark contribution to the D 
meson asymmetry, A corr , is small, less than 5%, compared to the measured asymmetries. The gluon 
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polarisation at NLO accuracy, obtained as weighted average over all bins, is: 

A^\ NL ° 

— > = -0.13 ± 0.15 (stat.) ±0.15 (syst.) (15) 
g I 

It is determined in the interval 0.12 < x < 0.33 with a weighted (x) rj 0.20, at the scale (/J. 2 ) « 13 
(GeV/c) 2 . The gluon momentum fraction x is taken from the simulations. 

For a given experimental acceptance for open charm tagging, the average value of x depends on the order 
of the QCD calculations used in the analysis. At the same time the energy in the photon-gluon centre- 
of-mass system required to produce a D° meson is higher with parton shower simulation as compared 
to the case of LO where parton shower is not simulated. Therefore a value of (x) at which the gluon 



polarisation is determined at NLO (jc) « 0.20, is higher than (x) 0.11, see for example Fig. 12 




Figure 12: Distributions of the gluon momentum fractions x for the simulated Djj^ events at the LO accu- 
racy (marked 'AROMA') and at LO with parton shower (AROMA + PS). Note different normalisations 
of the samples. 



A systematic uncertainty of the result in Eq. ((T5|) is estimated as follows, see Table 13 Contributions 



from P^, P u f, s/{s + b) and D are the same as discussed in Sec. 4.2 i.e. 5%, 5%, 2%, 7% and 1.6% 



respectively. Contributions from false asymmetries and the assumption given in Eq. ([9]), the same as in the 



statistically optimised method at LO, are discussed in Sec 5.2 Concerning the partonic asymmetry <zll, 
the following contributions resulting from the NLO determination of aix were studied: its dependence on 
scale, the charm mass value and the Monte Carlo mode (with or without PS). The renormalisation scale 
(here chosen equal to the factorisation one, Ref. [41 J) varied from m c to 3m c . This changes the gluon 
polarisation at most by a factor of two, leading to a conservative contribution to the systematic uncertainty 
of 0.1. The variation of the charm quark mass between 1.3 and 1.6 GeV/c 2 results in a contribution of 
0.05. The systematic uncertainty contribution from the simulation method was estimated using AROMA 
with and without parton showers. Note, that measured D meson spectra are described reasonably well by 
AROMA both with and without parton showers, and the difference in «ll between both cases is mainly 
due to the different phase space available for 'unobserved' partons. Such comparison allows us to give 
a conservative estimate of uncertainty due to the simulation method equal to 0.04. The total systematic 
uncertainty of (Ag/g) NL0 is obtained by adding all the contributions in quadrature and amounts to 0.15. 
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Table 13: Contributions to the systematic uncertainty of {Ag/g} 



Source 


S ((Ag/g)) 


Source 


sm/g)) 


Beam polarisation 


0.006 


s/(s + b) 


0.009 


Target polarisation P t 


0.006 


aLL 


0.119 


Dilution factor / 


0.003 


False asymmetry 


0.080 


Assumption, Eq. (|9j) 


0.025 


Depolarisation factor D 


0.002 


Total uncertainty 0. 146 



The result on (Ag/g) NLO , Eq. ( 15 ), was included in NLO QCD fits of polarised parton distributions, see 
Appendix for details. The fitted distributions of Ag(x)/g(x), evolved to Q 2 = 13 (GeV/c) 2 , are shown in 



Fig. 13 together with error bands corresponding to the statistical errors as derived from the error matrix 
of the fitted parameters. The present NLO open charm result agrees within 0.5 standard deviation with 
the fitted COMPASS curve for Ag(x) < and within 2 a with the one for Ag(x) > 0. It significantly 
influences the Ag(x) > fit, reducing the value of AG from 0.39 ± 0.07 (stat.) to 0.22 ±0.08 (stat.) at Q 2 
= 3 (GeV/c) 2 . 

The results of two other global fits, DSSV |42) and LSS 02, which employ both DIS and SIDIS asym- 
metries are also shown in the same Figure. In the DSSV fit, Ag(x) changes sign at x « 0.1 and is about 
1.5 a above the COMPASS open charm value. In the case of LSS, two solutions, with positive and with 
sign-changing Ag(x) are quoted. The LSS fit cannot distinguish between the positive and sign-changing 
Ag(x) functions. Both solutions give a positive Ag(x) at the (x, Q 2 ) of the present measurement, about 2 
and 1.3 a above the measured value. 

0.£ 



.O) 0.6 h 

< 0.4H 



Q 2 = 13(GeV/c) 2 



■k COMPASS open charm 

COMPASS (Ag>0) 

COMPASS (Ag<0) 

DSSV 

LSS(Ag>0) 

LSS (Ag sign changing) 



10* 



10"' 



Figure 13: The present NLO measurement of the gluon polarisation (Ag(x) /g(x)) at {ji 2 ) = 13 (GeV/c) 2 , 
compared to the NLO QCD fits of COMPASS with Ag(x) > (continuous line) and Ag(x) < (long- 
dashed) with their respective error bands, of LSS 11121 (dashed and dotted curves, respectively with 
Ag > and Ag changing sign) and of DSSV l42l (dashed-dotted curve), all at the same value of Q 2 = 13 
(GeV/c) 2 . The measurement error and the error bands are statistical only; the horizontal bar marks the 
range of x in which (Ag(x)/g(x)) is determined. 



6 Conclusions 

We have presented new results on the gluon polarisation in the nucleon {Ag/g) and the virtual photon- 
nucleon asymmetries A rN obtained from charm production tagged by D meson decays in 160 GeV/c 
polarised muon scattering off polarised proton and deuteron targets. The results are based on a data 
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sample collected between 2002 and 2007 and supersede the previously published ones IU91 as they are 
based on an enlarged data sample and an improved analysis method. The gluon polarisation is determined 
with open charm production dominated by the photon-gluon fusion mechanism followed by a spin- 
independent charm quark fragmentation into D mesons. This analysis neglects any contributions from 
the intrinsic charm; a contribution of resolved photon interactions was found negligible. 

Only one charmed meson is required in every event. This meson is selected through its decay in one 
of the following channels: D* (2010)+ -)• D°7r s | ow -)• (K7r/K7r7r°/K7r^7r)7r s | ow or D° ^Kk . The decays 
are selected using the invariant mass distributions of identified kaons and pions. A neural network is 
used to distinguish signal from background events in the data. The asymmetries A rN are extracted from 
these open charm events in bins of D° transverse momentum and laboratory energy. The average gluon 
polarisation obtained from these asymmetries at LO QCD accuracy amounts to (Ag/g) LO = — 0.10± 
0.22 (stat.) ±0.09 (syst.). This result is confirmed and statistically improved by employing a statistically 
optimised method of extracting the gluon polarisation: 




= -0.06±0.21 (stat.) ±0.08 (syst.). 



Both results are obtained in the range 0.06 < x < 0.22 of the gluon momentum fraction with (x) w 0.11 
and a scale (jU 2 ) « 13 (GeV/c) 2 . 

For the first time in this analysis, next-to-leading order QCD calculations for the determination of the 
gluon polarisation are employed, sing asymmetries A rN the gluon polarisation is obtained as: 




= -0.13 ±0.15 (stat.) ±0.15 (syst.). 



In this case the range of x, 0.12 < x < 0.33, leads to a higher average value, (x) rj 0.20, while the scale 
is approximately the same, ({J, 2 ) « 13 (GeV/c) 2 . 

The present measurement at LO QCD accuracy of the gluon polarisation in the nucleon, together with 
other measurements of SMC, COMPASS and HERMES, all situated around x ~ 0.1, points towards a 
small gluon polarisation at that value of x. This is a hint for a small value of the first moment, AG, of the 
gluon helicity distribution, although it in principle does not exclude a large value of AG. 

The {Ag/g}^ result was included in NLO QCD fits of polarised parton distributions. It significantly 
influences a fit in which Ag(x) > was assumed, reducing the value of AG from 0.39 ± 0.07 (stat.) to 
0.24 ± 0.09 (stat.) at Q 2 = 3 (GeV/c) 2 , after it is included. 



Appendix 

In this appendix, new NLO QCD fits of polarised parton distributions, including the (Ag/g) N result of 
Eq. ( [T5] > are presented. In our previous fit [3 ] the gluon helicity distribution is parametrised at a reference 
Q 2 of 3 (GeV/c) 2 in the form 

Ag(x) = ^ > . (16) 

x a H\-xf*&x 



Here r\ g is the integral of Ag{x), r\ g = AG. The same parameterisations, Eq. ( 16 1, are used for the singlet, 
non-singlet quark and gluon helicity distributions except for the singlet quark in the fit with r\ g > where 
a factor ( 1 + yx) is added in order to allow a change of sign. The high x parameter of the gluon helicity 
distribution is fixed to /3 ? = 10 in this r\ g > fit so that the total number of free parameters remained 
equal to 10. Two fits, in NLO QCD approximation and using all inclusive data with Q 2 > 1 (GeV/c)" are 



performed: one with Tj g > 0, the other with r\ g < 0. Both of them gave a comparable % 2 probability. 
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In the new fit all the data used in Ref. O are employed with an addition of the 15 COMPASS values 
of A p x published later B3l . As in the previous fit only statistical errors are considered. The reference 



Q 2 is kept at 3 (GeV/c) 2 and the same parameterisations, Eq. ( |16| >, are used. The total number of free 
parameters is also equal to 10. 

The new open charm result is not attached to a precise value of x and thus its contribution is taken into 
account by the average 



l /-0.33 

<*g> 



(0.33- 0.12) Jo.n 



^W = 13) 



dx (17) 



which is re-evaluated during the fit for any modification of one of the gluon or singlet quark parameters. 
The obtained value of (R g ) is compared to the open charm result voc = —0.13 with the statistical error 
0bc = 0.15 and the % 2 of the fit is incremented by ((Rg) — voc) 2 /°oc- 



The unpolarised gluon distribution g(x, Q 2 ) in Eq. (17i is taken from the MRST04 parameterisation 
ll44l . It was also used in Ref. [3 ]. In contrast to previous parameterisations of the same group, MRST04 
predicts a slower decrease of the gluon distribution at high x, (1 — x)$ with /3 ~ 3-4. For this reason 
the choice of /3 g = 10 for the fit with T] g > in Ref. Q leads to a strongly peaked distribution of Ag/g 
which in turn generates a dip in the fitted distribution of g^x) around x = 0.25 for low values of Q 2 
and leads in some cases to very asymmetric errors due to the limits imposed by the positivity condition 
|A?( X )I < g(x). To avoid these unphysical features, /3 g is now fixed to 6 in the fit with Tj g > 0. 

The present open charm result has practically no effect on the fit for Ag(x) < 0, where r\ g = —0.34 ±0.12 
with- and without that measurement, while it reduces significantly the positive T] g , from r\ g = 0.39 ± 0.07 
(stat.) to T] g = 0.22 ±0.08 (stat.) at Q 2 = 3 (GeV/c) 2 , after it is included. Similarily the values of a g for 
the fits including the open charm point are a g = 1.31 ±0.47 (stat.) for Ag(x) > and a g = 0.26 ±0.48 
(stat.) for Ag(x) < 0, both at Q 2 = 3 (GeV/c) 2 . 
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